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ABSTRACT 

Using a stellar mass limited sample of ~ 46, 600 galaxies (M* > 10 91 Mq) at 0.5 < 2 < 2, we show 
that the stellar mass, rather than the environment, is the main parameter controlling quenching of 
star formation in galaxies with A/* > 10 10 Mq out to 2 = 2. On the other hand, the environmental 
quenching becomes efficient at 2 < 1 regardless of galaxy mass, and it serves as a main star formation 
quenching mechanism for lower mass galaxies. Our result is based on deep optical and near-infrared 
imaging data over 2800 arcmin , enabling us to negate cosmic variance and identify 46 galaxy cluster 
candidates with M ~ 10 14 Mq. From M* ~ 10 9,5 to 10 10 ' 5 Mq, the fraction of quiescent galaxies 
increases by a factor of ~ 10 over the entire redshift range, but the difference between cluster and 
field environment is negligible. Rapid evolution in the quiescent fraction is seen from z = 2 to z = 1.3 
for massive galaxies suggesting a build-up of massive quiescent galaxies at z > 1.3. For galaxies with 
M* < 10 10 Mq at z < 1.0, the quiescent fraction is found to be as much as a factor of 2 larger in 
clusters than in field, showing the importance of environmental quenching in low mass galaxies at low 
redshift. Most high mass galaxies are already quenched at z > 1, therefore environmental quenching 
does not play a significant role for them, although the environmental quenching efficiency is nearly 
identical between high and low mass galaxies. 

Subject headings: galaxies: clusters: general — galaxies: high redshift — galaxies: evolution 
galaxies: stellar content — galaxies: star formation 


1. INTRODUCTION 

Galaxies define basic separate, independent entities in 
the universe and are the building blocks of her. There¬ 
fore, the study of their evolution throughout the history 
of the universe is the very core in our understanding of 
the universe. 

The A cold dark matter (CDM) cosmological models 
say that galaxies form through the gravitational collapse 
of baryonic matters — mostly in the form of gas — inside 
dark matter (DM) halos. As gas collects within DM ha¬ 
los, it begins to cool and form stars, which gives galaxies 
the very glittering look as we observe them. An inter¬ 
esting question regarding star-forming (SF) activity of 
galaxies is how this activity evolves and what governs 
this activity. Unlike the collapse of DM, which is mainly 
governed by gravity, SF activity is governed or affected 
by various baryonic physical processes - including cool¬ 
ing of gas as well as feedbacks from stars, supernovae 
(SNe), or active galacitc nuclei (AGN) — thus making a 
detailed understanding of the SF activity complicated. 

Regarding the evolution of SF activity of galaxies, 
we know that galaxies in the local universe have much 
lower levels of SF activity than the ones in the past. 
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First, the global SF hi story, revealed from various 
galax y surveys dHopkins fe B cacom 2006; Cuccia,ti et al.1 
2012t iBehroozi et al.1 20131 iBurgarella et al.1 120131 


Magnelli et al.l 120131 : iMadau fc Dickinsonl 1201411 shows 


significant drop from redshift, x ~ 1 to 2 ~ 0 — with 
a possible broad peak around 2 ~ 1 — 3. While this 
shows the behavior of the collective SF activity in the 
universe, we can still speculate from this that the SF 
activity in individual galaxies must also decrease signif¬ 
icantly. Second, star-formation rates (SFRs) of normal 
SF galaxies at 2 ~ 0 are significantly lower than SFRs 
of hig her redshift (2 >1) SF galaxies with similar mass 
(e.g.. iDaddi et al.l 120071 ISobral et aPl2014l l. This sug¬ 
gests an average drop in SF activities of individual galax¬ 
ies with time toward 2 ~ 0. Lastly, a large fraction of 
local galaxies — especially massive ones — are red, quies¬ 
cent, and form a tight red sequence (e.g . iStrateva et all 
120011 : iBlanton et al.l 120031 iBalogh et al.i 12004( 1 . There- 
fore, we need to understand the physical mechanism(s) 
that can lower and eventually stop the SF activity of 
galaxies to explain this change of average SF properties 
of galaxies from high redshift to the local universe. 

At low redshift, it has been well known that SFR, color 
and morphology of galaxies, show strong dependence on 
their mass (e .g. iKau ffmann et al.1 120031 I Jimenez et al.1 
12003 iBaldrv et. alfeooolh in such a way that less massive 
galaxies show more SF activities than massive ones. This 
clear mass-dependence seems to persist at higher red¬ 
shift. By analysing the red-sequence luminosity function 
of a la rge sam p le (~ 500) of galaxy clusters at 2 < 0.95, 
iGilbank et al. ' (l2008t l found a dearth of red galaxies at 
faint luminos ity at higher redshift. R ecentl y, using the 
UltraVISTA dMcCracken et al.l I2012H data, lllbert et al.l 







































2 


Lee et al. 


(120131) have shown the stellar-mass dependence in SF 
quenching at z < 1 in a sense that stellar mass density 
of massive (> 10 112 A/q) quiescent galaxies remains un¬ 
changed while less massive galaxies keep being quenched 
i.e., massive ones already quenched before z ~ 1, 
earlier than less massive ones. The origin of this mass- 
dependent SF cessation or quenching is still uncle ar. It 
could be the effect of feedbac k from AGN ( Hopkin s et al.l 
120061: [Somerville et al.ll2008h . which is believed to selec¬ 
tively act for more massive galaxies. Another compelling 
expla nation is the heating of ac creted cold gas in massive 
halos (jBirnboim fc Dekel120031 ). 

On the other hand, the SF activity — or related prop¬ 
erties, such as color or morphology — shows environ¬ 
mental dependence as well. Locally, galaxies show dis- 
tinct SF properties and colors in different environment 
dLewis et al. 1120021: iKauffmann et al1l20f)4l iBlanton et al.1 

120051) . Naturally arising questions are how and when 
these environmental trends have been developed. By 
analysing galaxy colors o ut to z ~ 1.5 from VIMOS-VLT 
Deep Survey (VVDS), iCucciati et al.l (|2006l) found that 
the color-density relation progressively weakens with in¬ 
creasing redshift and pos sibly reverses at 1.2 < z < 1.5. 
IGriitzbauch et al.1 (|2011h . based on th e Palomar Ob¬ 
serva tory Wide-Field Infrared (POWIR; iConselice et al.1 
120081) survey data in the DEEP2 field, also argued that 
there is no or weak correlation between galaxies’ local 
number density and the color or the blue galaxy fraction 
for 0.4 < z < 1 galaxies. I n the case o f the relation be¬ 
tween SFR and density, ICooper et al.l (I2008D found the 
reversal of the relation between SFR and overdensity at 
2 ~ 1 using galaxy samples from the DEE P2 G alaxy 
Redshift Survey (GRS). With GOODS data. lElbaz et al.l 
: 2007 ) also argued that the relation between SFR and 
density was reversed at z ~ 1, showing a sharp contrast 
to the local trend. 

On the contrary, ICooper et al.1 ([2006P argued that the 
mean galaxy environment shows a strong dependence on 
galaxy color at z ~ 1, similarly with what is locally 
found. [Patel et all (2 009f ). analysing z ~ 0.8 galaxies, 
found that the total SFR local-density relation still per¬ 
sists at this redshift. This environmental dependence 
of galaxy properties has been shown to persist even at 
higher redshift {z < 2) bv lQuadri et al.l (12012I) . through 
the analysis of quiescent galaxy fraction of UKIDSS/UDS 
galaxies, even though the environmental trend of qui¬ 
escent fraction becomes weaker with redshift. Com¬ 
bining data L Jrom_seyeral_surveys — UltraVISTA, 3D- 
HST dBrammer et al.l I2012D . Cosmic Assembly Near- 
infrared Deep Extragalactic Legacy Survey (CANDELS; 
Grogin et al.1 120111 : iKoekemoer et al.l I2011D . and SDSS, 
Tal et al.1 (120141) have studied the evolution of quiescent 
galaxy fraction in groups over a wide range of redshift 
up to z < 2.5. Their finding is that quiescent fraction of 
group satellite galaxies is similar with background (i.e. 
field) galaxies at z ~ 2, while it increases faster than 
field galaxies at lower redshift down to z ~ 0. So, there 
seems to exist a certain degree of disagreement on how 
and how much the environmental trend changes as we go 
up to high redshift, even though there is a consensus on 
the fact that the environmental depen d enc e weakens as 
redshift increases (e.g., iScoville et ahll2013D . Therefore, 
we have not yet reached to a firm conclusion about the 
rate or amount of evolution in the effects of environment 


on the formation and evolution of the quiescent galaxy 
population. 

In studying the role of environment in formation and 
evolution of quiescent galaxies, galaxy cluster offers a 
very unique laboratory. They define the densest galaxy 
environment, thus are expected to follow the highest den¬ 
sity peaks in the universe. Besides the fact that denser 
structures collapse earlier in the ACDM paradigm, galax¬ 
ies in these dense environments would also be subject to 
sever al physical processes — such as ram-p ressure strip¬ 
ping (IGunn fe Gotti 1 19721) . strangulation ([Larsoi^et^ 
1980; Balogh et alJl200(3)7 and harassment (IMoore et al 


19981 ) — acting exclusively in group- or cluster-like envi¬ 


ronments. In this regard, studying the stellar population 
properties of galaxies in galaxy clusters at various red- 
shifts is a useful and crucial test bed in the investigation 
of the environmental effects on the formation and evolu¬ 
tion of quiescent galaxy population. 

While there is a sharp contrast in galaxy properties 
between galaxies in clusters and in field locally, as the 
cluster environment is d ominated by the red, quies¬ 
cent galaxies (e.g.. IDressletill980l) . th is picture seems to 
change as we go out to high redshift (IBut cher fc Oemlerl 
11984 IDressler et al.l [1997 : IPoggianti et al l 120011) . Conn 
bining their field sa m ple with the cluste r samp les from 
iHolden et al.1 (I2007D . van der Wei et al.l (120071) showed 
that the difference in the red galaxy fraction between 
the cluster- and the field environments increases with de¬ 
creasing redshift i.e., the red galaxy fraction increases 
toward low redshift in the cluster environment while it 
remains nearly constant in the field environment. At 
higher redshift, 2 > 1.5, the story seems to become more 
complicated: Studies suggest that the galaxy SF prop¬ 
erties (the SFR or the color) in the cluster environment 
at z > 1.5 can be different from those of cluster galaxies 
at lower redshift and that there can be cluster-to-cluster 
variation in t he properties of these high-redshift proto¬ 
clusters. iGobat et al.l (|201lD analysed colors of red galax¬ 
ies in a proto-cluster at z ~ 2.1 originally identified by 
their red Spitzer /IRAC colors, and find that these red 
galaxies have old stellar populations. On the other hand, 
iZirm et ahl 1 2008 ) found a well developed red-sequence 
in a z ~ 2.2 proto-cluster around a radio galaxy and 
suggested that many of these red galaxies are p o rbab ly 
dusty star-forming galaxies. iFassbender et al.1 (120111) 
also reported the existence of the actively star-forming 
galaxies in the X-ray selected proto-cluster at z ~ 1.6. 
iTran et all (2 01'ct ) studied the color and SF properties of 
the member galaxies of Spitzer -selected z ~ 1.6 cluster 
([Papovich et al.l[2C)10D . They also found that the galaxies 
in this high redshift cluster, unlike the ones in lower red¬ 
shift, are more dominated by blue galaxies. This, com¬ 
bined with the higher star-forming galaxy fraction within 
this cluster than in the lower redshift (z ~ 0.3), indi¬ 
cates that many cluster galaxies ar e sti ll activel y form¬ 
ing stars at z ~ 1.6. The results of iHilton et al.l (l2010i) . 
based on the X-ray and mid-IR (MIR) observation, also 
show the strong SF activity of the clust er member galax¬ 
ies at z ~ 1.5. iStrazzullo et al.l (2013 ) have found that 
there are massive star-forming galaxies along with pas¬ 
sive ones even in the core regio n of the IRAC-sel ected 
z ~ 2.1 proto-cluster studied by (jGobat et al.ll201l1) . All 
together, these observational results point to the redshift 
range of z > 1.5 as the active formation era of the galaxy 
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clusters as well as the epoch when many of cluster galax¬ 
ies are still forming stars actively. 

Summarizing, both the mass and environment of galax¬ 
ies seem to affect the quenching of SF activity of galaxies 
and the formation of quiescent galaxies. Also, the envi¬ 
ronmental dependence of SF activity of galaxies seems 
change with redshift. Therefore, more investigation 
about the timing (i.e., the redshift dependence) as well 
as the relative contribution from these distinct mech¬ 
anism^— mass v ersus environment — is still needed. 
iPeng et HI (|2010H have p rovided interestin g results in 
this r egard. Using SP SS ()York et all 12001)1) and zCOS- 
MOS ([Lilly et al.ll2007f) data, they analysed the effects of 
mass as well as environment on the evolution and star- 
formation activity of galaxies at redshifts up to z ~ 1. 
They found that “mass quenching” is more dominant for 
high stellar-mass (> 10 10 M 0 ) galaxies, while lower mass 
galaxies are more affected by “environmental quench¬ 
ing” . They also suggested that the combination of these 
two quenching mechanism can explain the evolution of 
stellar mass function of star-forming as well as quiescent 
galaxies. Analysing z ~ 1 Ha emit ters f rom the High¬ 
er Emission Line Survey (HiZELS), [Sobral et al.l ( 1201 1) 
have shown that the median SFR depends on galaxies’ 
environment for galaxies with their stellar mass lower 
than 10 10 ' 6 Mq, while there is no environmental depen¬ 
dence for more massive galaxies. 

In this work, we extend this kind of investigation to 
higher redshift (up to z ~ 2), examining the effects of 
these two suggested drivers - i.e. mass and environment 
in the formation and the evolution of quiescent galaxy 
population. It has been know n that a tight red sequence 
alrea dy exists at z ~ 1 (e.g. Ilm et al.ll2002t iBell et all 
200 1;. Therefore, investigation at higher (z > 1) red¬ 
shift is crucial in catching the formation of these quies¬ 
cent galaxies to place a meaningful constraints on their 
evolution. Our work is based on the deep near-infrared 
(NIR) data from the United Kingdom Infrared Telescope 
(UKIRT) Infrared Deep Sky Survey (UKIDSS) Ultra 
Deep Survey (UDS) as well as deep optical data from 
the Subaru/XMM-Newton Deep Survey (SXDS), com¬ 
bined with rnid-IR (MIR) data from the Spitzer Space 
Telescope. Deep UDS NIR data are crucial in robust es¬ 
timation of SF or stellar population properties of galax¬ 
ies, breaking the degeneracy between old and dusty star¬ 
forming populations, thus in a study of quiescent galaxy 
population evolution at redshift as high as z ~ 2. Our 
analysis is based on (1) an un-biased cluster detection 
(unlike, for example, the red-sequence method), (2) a 
relatively long redsliift-baseline (0.5 < z < 2.0) and a 
large spatial area (~ 665 Mpc 2 at z ~ 1), and (3) robust 
estimation of SFR and stellar mass from SED-fitting. As 
demonstrated below, deep NIR data in the UDS field en¬ 
ables us the robust estimation of the stellar population 
properties of the galaxies up to z ~ 2. 

We present the data set used in this work as well as 
the sample selection procedure in Section 2. In Section 
3, we explain how the photometric redshifts as well as 
the stellar population properties of the UDS galaxies are 
estimated. We explain the identifying procedure of the 
high-redshift clusters and provide the properties of the 
selected cluster candidates in Section 4. We analyse the 
stellar population properties and the evolutionary trends 


of the quiescent galaxies, focusing on the drivers of this 
population in Section 5, and we summarize our results in 
Section 6. We adopt the standard flat ACDM cosmology, 
with (U m , Ua) = (0.3,0.7), and H 0 = 70 km s -1 M pc~_. 
All magnitudes are given in AB magnitude system ( Okel 
11974) except when otherwise mentioned. 

2. DATA AND SAMPLE 

The UDS (Almaini et al. in prep.), which is one o f 
the five surveys in the UKIDSS (|Lawrence et all [2 00 71) . 
provides a very deep NIR imaging dataset over an area 
of ~ 0.77 degree 2 located on the SXDS. The survey 
was carried out with the Wide Field Camera (WFCAM; 
iCasali et al.1 l2007h on the UKIRT in three NIR broad¬ 
bands with the 5-er limits of J = 24.3, H = 23.3, and 
Ks = 23.0 (in Vega, DR10). 

The SXDS provides the deep optical data from the 
SUPRIMECAM on the Subaru telescope from B-band 
through z'-band. The 3-cr depths of the Subaru data 
are B = 28.4, V = 27.8, R = 27.7, i' = 27.7, and 
z' = 26. 7. We use the release d DR1 (z'-band detected) 
catalog (|Furusawa et al.l|2008f) . which is obtained from 
the SXDS DR1 release pag^l, for this work. The UDS 
region has also been observed by the Spitzer/TRAC as 
the SpUDS Spitzer Legacy Survey (PLDunlop), reaching 
~ 24 magnitude (Channels 1 and 2). We obtain the cata¬ 
log at the Spitzer Legacy Survey archiviQ, and aperture- 
corrected total magnitudes are used for Spitzer/TRAC 
data. We also use the publicly available spectroscopic 
redshift data for about 4000 objects (Simpson et al. in 
prep, Akiyama et al. in prep, and Smail et al. 2008) 

We perform photometry on the UDS DRIP i mages 
( J,H,Ks ) using SExtractor dBertin fc Arnoutsl Il996f) 
software. We use 9x9 convolution mask of a Gaus¬ 
sian PSF with FWHM of 5 pixels, and set the detection 
minimum area as 6 pixels. The detection and photome¬ 
try are done in each band separately, and then the match 
between different bands is done with a matching radius of 
1” and with the Tvs-band image as the reference image. 
About 0.8% objects (both in J- and 77-band catalogs) 
have multiple matches, and the nearest one is selected as 
the best match. To construct the panchromatic spectral 
energy distributions (SEDs) from the optical to the MIR 
in the observed frame, we match, first, the Subaru op¬ 
tical and UDS NIR catalogs, and then with the SpUDS 
IRAC catalog, with a matching radius of 1" and the posi¬ 
tion in A's-band catalog is used as the reference position. 
For optical and NIR data, we use the auto magnitudes, 
and the aperture-corrected total magnitudes are used for 
IRAC data. 

We check the validity of using auto magnitude 
( Mag au t 0 ) in crowded region, such as clusters, by com¬ 
paring auto magnitude values with aperture (circular di¬ 
ameter = l'/5) magnitude (Mag aper ) in held and cluster 
regions. If Mag au to values are biased in cluster regions, 
then we expect that Mag au t 0 to deviate from Mag aper 
differently in cluster and held. At several magnitude 
bins (Ks = 19.9, 20.9, 21.9, 22.9), we compare the distri¬ 
butions of Mag apel — Mag au to between cluster and held 
region through the Kolmogorov-Smirnov test (K-S test). 
At each magnitude bin, the maximum difference, D , was 

6 http://soaps.nao.ac.jp/SXDS/Public/DRl/ 

7 http://irsa.ipac.caltech.edu/data/SPITZER/SpUDS/ 
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0.073, 0.086, 0.071, 0.058, and the probability of null hy¬ 
pothesis, p, was 0.93, 0.62, 0.64, 0.94, showing that there 
is no systematic bias in auto-magnitude measurements 
that is specific in cluster region. 

In our analysis, we use objects which are detected in 
both of optical (Subaru) and NIR (UKIRT) catalogs as 
our sample, but regardless of IRAC detection. There are 
about 8.4% of Jvs-band detected objects not detected 
in z' band. The fraction is lower as 5.9% if we only 
consider the sources with Ks < 24.9. These z'-band 
undetected sources have, on average, redder IR colors 
than the detected ones. The mean values of (J — Ks) 
and (J — 3. 6/7777) are 1.0 and 1.8, respectively, while the 
corresponding values are 0.6 and 1.0 for z'-band detected 
objects. And, about 25% of Ks- and z'-bands detected 
objects are matched with IRAC sources. 

From this ‘Subaru+UKIRT+Spitzer’ multi-band pho¬ 
tometric catalog, we first cull out stars either based on 
the spectroscopic classification or broad-band color (BzK 
or J-K ). Specifically, from ~ 200,000 total objects, we 
exclude: (1) ~ 550 spectroscopically classified stars, (2) 
~ 22,000 stellar candidates based on their location in 
the ( B-z')-{z'-Ks ) color-color plane, with the color cri¬ 
terion of (z 1 — Ks) < 0.3 x ( B — z') — 0.4, s imilarly with 
iDaddi et all (1200411 or lHanami et al.l (120121) . and (3) fur¬ 
ther ~ 500 candidates with J — Ks < 0 and Ks < 19.0, 
based on the fact that galaxies and stars form clearly 
separate branches in this J — Ks versus Ks plane. 

We need to cull out active galactic nuclei (AGN) from 
our sample because the measurements of photometric 
redshift and stellar population properties, which are cru¬ 
cial in our study, would be affected by the AGN activ¬ 
ity. To remove AGN candidates from this star-removed, 
multi-band catalog, we use three steps: First, we ex¬ 
clude 77 objects which are classified as quasars in the 
spectroscopic catalogs. Next, we match our multiband 
catalog wit h the XMM-Newton point source catalog 
(jUeda et al.l 120081) to remove any matched object. At 
this step, about 220 objects are matched and excluded. 
Finally, for the objects with the matched IRAC photom¬ 
etry, we use the NIR & MIR color criteria suggested by 
iMessias et all ( 2012 ) to find out any remaining AGN can¬ 
didates. The applied criteria are ( K — 7714 . 5 ) > 0 and 
( 7774.5 — 7778. 0 ) > 0 , where, 7774.5 and 777 s .0 are the magni¬ 
tudes in IRAC channel 2 (4.5/7777) and channel 4 (8.0/ 7 m), 
respectively. About 7,000 objects satisfy these color cri¬ 
teria, thus are excluded from the sample. These color 
criteria, similarly with other IRAC color criteria, can¬ 
not distinguish between AGN and higli-redshift (z > 2) 
galaxies. Therefore, by applying these criteria, we might 
lose some non-active galaxies at z > 2 as well. How¬ 
ever, this does not affect our results, because our target 
redshift range does not exceed z ~ 2. 

After trimming out AGN candidates as well as stars as 
explained, we apply the A's-band magnitude cut (Ks < 
24.9, 5 cr detection limit) and the number of remain¬ 
ing galaxies is about 115,000 among which ~1,400 ob¬ 
jects have spectroscopic redshift information. For these 
~115,000 objects, we derive photometric redshifts as 
well as stellar population properties through SED-fitting 
(Section 3). After selecting galaxies within the redshift 
range of 0.45 < z < 2.1 and applying the stellar mass 
cut of log M^/Mq > 9.1 — which is the 75% stellar mass 
limit at the highest redshift bin (z ~ 2) of our study — 


our final sample size is 46,641. 
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Figure 1. (left) Comparison of spectroscopic and photometric 
redshifts for the subset of UDS (red) and GOODS-S (blue) galaxies 
with available spectroscopic redshift. (upper right) Photometric 
redshift discrepancy as a function of redshift. (lower right) Out¬ 
lier (|Az|/(l + z sp ec ) > 0.15) fraction as a function of redshift. 
Numbers above each bar show the number of outliers at each red¬ 
shift bin. 
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3. ESTIMATION OF PHOTOMETRIC REDSHIFTS AND 
STELLAR POPULATION PROPERTIES 

3.1. Photometric Redshift 

We derive photomet ric redshifts using the EAZY code 
(IBrammer et al.l I2008D . We take the AT-band magni¬ 
tude prior, which slightly improves the photometric red¬ 
shift accuracy. We fit ten broad-band photometric data 
points from the observed frame optical (Subaru) to MIR 
(Spitzer/ IRAC). During the fitting procedure, we re¬ 
quire that at least five colors to be measured - i.e. it 
is detected at least at six bands. Therefore, any galaxy 
which is not detected at more than five bands is not in¬ 
cluded in the following analysis. Only ~ 0.14% of objects 
are excluded for this reason, and these are mostly faint 
objects — AT > 24 for ~ 72% among these, and zJ > 26 
for ~ 52%. 

For a subset of galaxies for which spectroscopic 
redshifts are available and Ks < 24.9, we compare the 
spectroscopic and the photometric redshifts to estimate 
the accuracy of the photometric redshift (Figure |I] 
left panel). In UDS, there are not enough galaxies 
with z spec > 1-5 for a fair estimation of photometric 
accuracy. To test the photometric redshift accuracy 
at z > 1.5, we use the GOODS-S data, which carries 
similar coverage of broad-bands data with the UDS while 
having more abundant spectroscopic redshift samples at 
z > 1.5. These GOODS-S spectr osco pic d ata ar e from 
various spectroscop ic sur veys dLe Fevre et al.l 12004 
Szokolv et, al.l 120041 iMignoli et al.1120051: iCimatti et al l 


20081: Vanzella et all 20081: iPopesso et al.1 120091 : 

Balestra et al. I 120101 : iSilverman et al.1 I2010D . We 
measure signal-to-noise of UDS photometric data as a 
function of magnitude at each band, and add additional 
noises (randomly following a Gaussian distribution) to 
the GOODS photometry to simulate the depth of the 
UDS data. Then, photometric redshifts are derived 
using these noise-added photometry. Within the redshift 
range of 0.4 < z < 2.1, the redshift range of interest 
of this work, there are 3% of objects whose redshift 
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error, defined as |Az|/(l + z spec ), is greater than 0.15. 
Removing these objects with catastrophic redshift 
errors, the mean error is 0.028. 

The right panels of Figure [H show |A^|/(1 + z spec ) as 
well as the outlier fraction as functions of photometric 
redshift. We can see that our photometric redshift esti¬ 
mation is reliable within the redshift range of our study, 
0.5 < ^ < 2, but its uncertainty is larger at z > 1.5. The 
outlier fraction increases at high redshift (z > 1.5), but it 
is kept at < 0.2. The reason for the increase in the photo¬ 
metric redshift uncertainty at z > 1.5 is that the spectral 
break at ~ 4000 A one of the clearest redshift indica¬ 
tors — moves out of the optical regime at z > 1.5. When 
that happens, it becomes more difficult to sample the 
exact location of the rest-frame 4000 A break in the ob¬ 
served frame due to the gaps between NIR filters and the 
reduced number of available filters above 4000 A break. 
Also, the fraction of star forming galaxies are higher at 
2 > 1.5 in the spectroscopic sample since it is easier to 
spot emission lines than absorption lines for faint galaxy 
spectra. Weak 4000 A breaks and strong emission lines 
of star forming galaxies make it difficult to trace the 
continuum shape, leadi ng to a reduced accuracy in pho¬ 
tometric redshifts (e.g.. rD amen et al.||2009|: [Ilbert^t^d 


2008 ; IQuadri et all 20121 IHartlev et alJl20l4lYang et aT 


2014 iKang fc Iml 120151) . Not surprisingly, photometric 
redshift outliers tend to be star-forming galaxies, and 
some of them can move into the red sequence due to the 
large errors in photometric redshift. 

In Figure 2, we show the z spec distribution of galax¬ 
ies with several z p hot bins (z ~ 0.6,0.8,1.0,1.2,1.4 and 
1.6) with |As| < 0.1. Figure 2 shows that systematic 
offsets between photometric and spectroscopic redshifts 
are small (< 0.1 at some redshift bins). The effect of 
the small systematic shift in redshift is negligible in our 
analysis. For example, a systematic offset of A z ~ 0.1 
causes a shift in the stellar mass estimate by < 20%. 



0.5 1.0 1.5 2.0 0.5 1.0 1.5 2.0 0.5 1.0 1.5 2.0 



Figure 2. Distributions of spectroscopic redshift of galaxies se¬ 
lected by their photometric redshift at various redshift bins. For 
example, at redshift 0.6, this figure (upper left panel) shows the 
distribution of spectroscopic redshifts of galaxies whose photomet¬ 
ric redshifts are within 0.5 < z p h 0 t < 0.7. Systematic offsets are 
small, if there is any (j 0.1), between photometric and spectro¬ 
scopic redshifts. The fraction of redshift outlier increases at higher 
redshift while the number of ourliers remain similar. 


3.2. Stellar Population Properties from SED-fitting 

We perform the stellar population analysis through 
spectral energy distribution (SED) fitting methods for 


our sa mple, with t he same pro cedure as explained in de¬ 
tail in ILee et all (120091 1201Of) . For redshift, we adopt 
z apec if available, but z p hot otherwise. 

The stellar p o pulatio n synthesis models of 
iBruzual fc Charlotl (j2003l BC03, hereafter) with 

Padova 1994 stellar evolutionary tracks are used for our 
SED-fitting analysis. We assume the iChabrieJ (120031 ) 
initial mass function (IMF) with the lower and upper 
mass cuts of 0.1 and 100 M 0 . 

As for the star-form atio n hist ories (SFHs), we use the 
delayed SFHs (ILee et al.l [20101 1. This form of SFH is 
shown to provide better estimation of stellar population 
properties than ot her forms, for example, exponentially 
declining ones (e.g.l Marast on et al.ll2010l : iPapovich et al.1 
l20HlPforr et~aITl2012fc ILee et al.ll2014ll . Stud ying 2 ~ 1 
galax i es from DEEP2 Red shift Survey ([Newman et al.1 
120 1 ,‘lt ) . iPacifici et al.1 (|2013H have also shown that these 
galaxies have SFHs similar with this form. 

It has the following functional form. 

$(t,r)cc4 e " i/T , (!) 

where, 'I'(t, r) is the instantaneous SFR for a given set 
of r and t. The parameter t is the time since the onset 
of the star-formation and r is the time-scale parameter 
which governs how fast (or slowly) the SFR reaches its 
peak value before starting to decline. We allow the value 
of r to vary within a very broad range from 0.1 Gyr to 
10.0 Gyr, and t within 0.1 Gyr < 2 < tjj(z), where tjj(z) 
is the age of the Universe at redshift 2 . 

After extracting ~ 4,000 spectral templates from 
BC03 with t hese allowed values of parameters, we apply 
the iGalzetti et al.l ( 2000 ) attenuation law to model the 
dust attenuation of SEDs due to the inter-stellar dust. 
We vary the amount of attenuation, parametrized via 
E(B — V ) within 0.0 < E(B — V) < 1.50 with the in¬ 
crement of 0.025. For the attenuation due to the line- 
of-sight neutral hydrogen in the inter-g alactic m edium 
(IGM), we apply the extinction law of iMadaul (119951 1. 
During the fitting, the redshift is fixed at the value of 
the spectroscopic redshift if available, otherwise at the 
photometric redshift. 

We apply the ‘grid-scanning’ meth od w hich has been 
used in several previous works dWiklind et al.l 120081 : 
ILee et all 12009112010112014H . As explained in these ref¬ 
erences, our method of skimming the entire parameter 
space is less subject to the effects of any prior. From 
this SED-fitting results, we derive for each galaxy, the 
stellar mass (M»), the SFR, the mean age of the stellar 
population, the amount of dust attenuation (E(B — V)), 
as well as the star-formation history — parametrized in 
terms of r and t. The SFR used in this work is defined 
as the SFR averaged over recen t 100 Myr, based on the 
reasoning outlined in ILee et al.l (1200911 . We also derive 
the rest-frame u — g color — will be denoted as (u — g )0 
from now on — for each galaxy from its best-fit BC03 
template. 

The broad wavelength coverage out to IR regime as 
well as the deep optical data from Subaru in the UDS 
field enables us to estimate the stellar population prop¬ 
erties of high-redshift galaxies with high confidence and 
distinguish dusty SF galaxies from old galaxies among 
red galaxies. Figure [3] shows the examples of the SED- 
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E(B-V) =0.125 
=300 Myr 



3.6 3.8 4.0 4.2 4.4 4.6 
log A [A] 


3.6 3.8 4.0 4.2 4.4 4.6 
log A [A] 
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log A [A] 


Figure 3. Examples of observed SEDs and the best-fit BC03 spec¬ 
trum for blue SF (left column), dusty SF (middle column), and 
old passive (right column) galaxies at 2 ~ 0.7, 1, 1.5. In each 
panel, red circles show the observed fluxes from BVRiz (Subaru) to 
JHKs (UKIRT) to MIR (Spitzer) upto IRAC Channel 2 (4.5 fim). 
Solid blue curve is the best-fit BC03 spectrum found through SED- 
fitting. We can see that blue SF galaxies have bluer colors than red, 
dusty SF galaxies throughout the shown wavelength range as well 
as that dusty SF galaxies and old quiescent ones can be discrim¬ 
inated reliably thanks to the deep and wide wavelength coverage 
of the data. Wavelength (rr-axis) is in logarithmic scale of A, and 
flux is in log-scale of 



Figure 4. Galaxy number counts (and error) at stellar mass bins 
(blue) and the same ones but with photometric redshift scatter 
effect included (red) at z = 1 (left) and z = 1.5 (right). Red solid 
and dotted lines in each panel show the mean and envelope (i.e. 
minimum and maximum) of Monte Carlo simulation (100 times) 
which include the redshift scatter. We can see that the effect of 
photometric redshift error is not significant in the estimation of the 
stellar mass function. 

shift (red) show no significant difference from the ones 
without scatter (blue) at both redshift bins. For indi¬ 
vidual galaxies, the mean and the standard deviation of 
is 0.02 ± 0.11 and 0.01 ± 0.09, at z ~ 1 
and ~ 1.5, respectively. Here, AM* is M* j mc — Af* jtrlie , 
and M* jtrlie is the measured stellar mass without redshift 
scatter. From this, we can ensure that the stellar mass 
is not significantly affected by the error in photometric 
redshift measurement. 


fitting results for a representative set of galaxies at var¬ 
ious redshifts. Also shown are the key SED-fit parame¬ 
ters. In this figure, we show the observed galaxy SEDs 
(red circles) as well as the best-fit BC03 spectra for typ¬ 
ical blue SF galaxies (panels in left column), dusty SF 
ones (panels in middle column) and old quiescent ones 
(panels in right column) at redshifts 2 ~ 0.7, 1, 1.5. 

Here, we define “quiescent galaxies” as galaxies with 
sSFR < l/[3t(z)\ yr -1 , where t(z) is the age of the uni¬ 
verse aHzMIhisajDwas used in several previous works 
(e.g.- lDamen et al1l2009j h and approaches the local value 
for sS FR cut at 2 ~ 0 (e.g. JGallazzi et aT1l2009l:iKo et, al.1 
120141 b Under a delayed SFH model (i.e., Equation (1)), 
l/[3t(2)] corresponds to ~ 1% of the peak SFR. While 
this sSFR cut select galaxies with low SF activity in com¬ 
parison to the average SF galaxies at that redshift, the 
selected galaxies can have SFRs that are not quite “qui¬ 
escent”. For example, at 2 ~ 2, M* ~ 10 11 Mq galaxies 
at the sSFR cut have ~ lOM 0 /yr — i.e., SFR of star- 
burst galaxies in the local universe. If we fit the SEDs 
of a 1 Gyr-old, passively evolving galaxies after a spon¬ 
taneous burst (a truly quiescent galaxy at high redshift), 
our procedure returns sSFR ~ 10” 11 yr -1 . We check 
the MIPS 24 /jm detection of our galaxies using SpUDS 
MIPS catalog. We find that small fraction (~ 5%) of 
red SF galaxies is detected, while no quiescent galaxy is 
detected at 24 /mr. 

We test how the error in photometric redshift mea¬ 
surement may affect the estimation of stellar masses of 
galaxies using Monte Carlo simulation. Specifically, first, 
we scatter the redshift of each galaxy assuming a Gaus¬ 
sian distribution of redshift error. Then we repeat SED- 
fitting by fixing the redshift at this scattered value. We 
repeat this procedure 100 times each galaxy at 2 ~ 1 
and 2 ~ 1.5, within a redshift range ±0.02. The re¬ 
sults of this simulation are given in Figure [2 where we 
show the number counts of galaxies at stellar mass bins. 
The resulting number distributions with scattered red- 


4. MASSIVE STRUCTURES OF GALAXIES 
4.1. Stellar Mass Limit 

Based on the stellar mass derived from the SED-fitting, 
we select galaxies with their stellar masses greater than 
10 9 - 1 M 0 (which is the 75% stellar mass limit at the 
highest redshift range, 2 ~ 2.0) within the redshift range 
of 0.45 < 2 < 2.1. 



log\ [A] 


Figure 5. Depths of the datasets used for our study, in compar¬ 
ison to different types of galaxies. Black horizontal bars show the 
flux limit of the UDS in 5 Subaru filters ( B , V, R , i' , z'), 3 UKIRT 
filters (J, H , Ks ), and 2 IRAC channels (3.6/xm, 4.6/xm). All the 
galaxy templates are normalized to have stellar mass of 10 91 Mq. 
Solid blue curve shows the spectrum of a young galaxy template 
with t = 200 Myr with no dust attenuation. Solid green curve is 
for the same galaxy template but with high dust attenuation value 
of E(B — V) = 0.4. Red curve shows an old, quiescent galaxy 
template with t = 2 Gyr — maximally possible age at z = 2 if 
we assume the formation redshift as Zf orrn = 5. This figure shows 
that we can detect with our data set galaxies at z ~ 2 with stellar 
mass as low as ~ 10 91 Mq. 

We have checked if the magnitude limit is deep enough 
to detect galaxies with stellar mass of 10 91 Mq out to 
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z = 2.0, using galaxy template spectra of various galaxy 
types extracted from the BC03 library. In Figure 0 we 
show these spectra — young (solid blue curve), dusty 
young (solid green curve), and old (solid red curve) - 
with stellar mass of 1O 91 M 0 and redsliifted to z = 2.0. 
We also show the limiting flux of the UDS data set in 5 
Subaru optical filters ( B , V, R, i', and z'), in 3 UKIRT 
filters (J, H 1 and Ks ) and in 2 IRAC channels (3.6/im 
and 4.5p,m) shown as black horizontal bars. We find 
that SEDs of all the galaxy types are above the detection 
limits in multiple bands, demonstrating that our dataset 
is nearly suitable for selecting galaxies with log (M*/M 0 ) 
> 9.1. 

Figure [6] shows the redsliift versus stellar mass of galax¬ 
ies in the Af-band limited sample. In this figure, the green 
vertical lines show the sample redshift range of this study 
and the green horizontal lines show the stellar mass lim¬ 
its. Blue curve shows the stellar masses at different red- 
shifts of star-forming galaxy templates from BC03 with 
formation redshift of Zf = 3. Yellow curve corresponds to 
quiescent galaxy templates with Zf = 5, which is a rea¬ 
sonable formation redshift for faint/low-mass galaxies, 
while red curve is for galaxy templates with extremely 
high formation redshift (zj = 10). As shown in this 
figure, within the given redshift range, we can detect 
log (M*/M 0 ) > 9.1 galaxies (corresponding lower green 
horizontal line) either quiescent (yellow curve) or star¬ 
forming (blue curve). It is still possible that we may miss 
some low mass (< 1O 95 M 0 ) galaxies with extremely old 
(zf ~ 10) population, even though it is reasonable to as¬ 
sume that there are only few, if any, galaxies with these 
extreme properties in reality. For this reason, we apply 
the stellar mass cut of 1O 9 5 M 0 (upper green horizontal 
line) when we deal with quantities which might be af¬ 
fected by this missing of low-mass old galaxies. Within 
this redshift range and above the stellar mass limit (log 
(M*/M 0 ) = 9.1), we have 46,641 galaxies. 

4.2. Method and Procedure 

Our method of finding galaxy cluster candidates is 
based on the ph otometric redshift and si milar with the 
method used in iKang fe iinl (|2009l I2015D . Following is 
the detailed explanation of the procedure we use to iden¬ 
tify the cluster candidates. First, we divide the projected 
sky area into retangular grids with a width of 12", which 
corresponds to ~ 96, 102, and 100 kpc at z = 1, 1.5, and 
2. Then, at a given redshift bin and at each grid point, 
we count the galaxies with log (M*/M 0 ) > 9.1 within a 
radius of r < 700 kpc (~ 1.5' at 2 ~ 1) from the grid 
point within the redshift range of the typical redshift er¬ 
ror in the UDS (i.e., A z = ±0.028 x (1 + 2)). We repeat 
this procedure at redshift bins with an increment of A z 
= 0 . 02 . 

Next, we find the mean ( N ) and the standard devia¬ 
tion (on) of the galaxy number counts ( N ) through the 
Gaussian fitting at each redshift bin. Then, we iden¬ 
tify the over-dense grid points with their galaxy number 
counts that exceed the 4 -0 level from the mean, i.e., 

N > IV ±4 x cr N . (2) 

If we see the spatial distribution of these over-dense 
grid points with extreme galaxy number counts (> 4cr), 
we find that some points are connected to other points 



0.5 1 1.5 2 

Redshift 



Figure 6. (Top) Redshift versus stellar masses derived from the 
SED-fitting. Green lines show the redshift range (vertical) as well 
as the stellar mass limits (upper horizontal: near 100% at z < 2; 
lower horizontal: near 100% at z < 1.5 and 75% at z < 2) of our 
sample. Solid curves shows the stellar mass at different redshift 
of quiescent (yellow and red) and SF (blue; Zf = 3) BC03 galaxy 
templates with K = 24.9. The yellow and red curves correspond to 
quiescent galaxies with reasonable (Zf = 5) and extreme (zj = 10) 
formation epoch, respectively. This figure illustrates that we can 
detect reasonably old galaxies with log (M»/Mq) = 9.5 up to z = 2, 
while we may miss some extremely old, low-mass galaxies, if any. 
(Bottom) Distribution of spectroscopic redshift. It shows peak 
near z ~ 0.6 and ~ 0.9 similarly as total (photometric redshift) 
sample as shown in top panel. The number of spectroscopic sample 
quickly drops at z > 1. 


making large structure while some points are relatively 
isolated. In our case, the spacing between the grid points 
(~ 100 kpc) is much smaller than the typical size of 
galaxy clusters (~ 1 Mpc). Therefore, a cluster of galax¬ 
ies would appear as a structure with connected over- 
dense grid points. For example, if 9 grid points are con¬ 
nected with each other with a square shape (i.e., 3x3) 
or if 10 grid points are connected in a rectangular shape 
(i.e., 2x5), this connected structure would represent 
an overdense structure with a radius of ~ 850 — 900 
kpc. Based on this reason, we exclude any over-dense 
grid points when the number of connected grid points is 
smaller than 10 at each redshift bin. Then, among the 
remaining connected structures, we only select the struc¬ 
tures as our final (proto-)cluster candidates only when 
these structures are identified at the same sky location 
(within the displacement, d < 1.5 Mpc) in at least three 
successive redshift bins. In other words, if the connected 
structure is found only in one or two redshift bin(s), we 
discard this structure. This criterion has been applied to 
minimize any false detection of cluster candidate. By ap¬ 
plying these two, rather conservative criteria in selecting 
cluster candidates, we try to make our selection of cluster 
sample as robust as possible. By comparing with a m ock 
catal og from the GALFORM simulation IlMerson et al.l 
HoH), we find that our cluster finding method is nearly 
complete at selecting Mh > 1O 13 - 8 M 0 clusters at z ~ 0.9 
and z ~ 1.5. A detailed description of this metho d ap- 
plied on simulation data can be found in IKang fc Iml 
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Redshift 


Figure 7. The number of identified cluster candidates within the 
redshift range of 0.5 < -z < 2.0. 

4.3. Galaxy Cluster Candidates 

Following the procedure as explained in the previous 
section, we find 46 cluster candidates within the red¬ 
shift range of 0.5 < z < 2.0. Among these, 19 clus- 
ters are the ones already identified by other authors 

(Ivan Breukelen et al.ll200dl2007t IFinoguenov et al.ll2C)10b 

iPanovich et al.ll2010fl in the same UDS-field, and 27 clus¬ 

ters are newly identified in this work. Figure [7] shows the 
number of the cluster candidates at each redshift bin. 

In Table 1, we present the basic properties of these clus¬ 
ter candidates, including the number of member galax¬ 
ies (Column 4), total stellar mass of member galaxies 
(Column 5), and the overdensity measure, defined as 
(N — N)/un (Column 6). 
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Figure 8. Total number of member galaxies (left panel) and 
the maximum overdensity factor (right panel) of our galaxy clus- 
ter candidates. Red ci rcles are 13 galaxy clusters identified by 
IFinoguenov et afl 1120101 ) . 


After finding the cluster candidates, we define the spa¬ 
tial center and the redshift of the candidate clusters 
through the following procedure. At each redshift bin, 
the mean position, weighted by the galaxy count (N), 
is found for each connected over-dense structure. Next, 
the center and the redshift of each cluster candidate is as¬ 
signed in a similar manner, i.e, as the N -weighted mean 
center and redshift, for each of the connected redshift 
bins. Then, we find the galaxies within the radius, r < 1 
Mpc — which corresponds to, or is slightly smaller than, 
a typical radius of a galaxy cluster — from the center 
and within the redshift interval, |Az| < 0.028 x (1 + z) 
- which is the typical redshift error of our data. Fi¬ 
nally, we re-calculate the center and the redshift of each 
cluster candidate as the stellar-mass weighted mean cen¬ 
ter and redshift of the member galaxies belonging to the 
overdense area. These center and redshift are the values 


listed in Columns (1), (2) and (3) in Table 1. The num¬ 
ber of member galaxies varies within the range between 
~30 and ~90 (the left panel of Figure[8]). The right panel 
of the same figure shows the over-density measure, de¬ 
fined as (N — N)/<jn , and shown as <jod (y-axis) in the 
figure. Most clusters have ood values between 4 and 8, 
and 4 clusters at redshift z < 1 have higher god- 
Out of 46 clusters, 13 clusters are also found by 
IFinoguenov et al.l (120101 1 based on the A'-ray observation, 
and they provide the total mass ( M 200 ) measured from 
the X-ray data, and the corresponding r 2 oo- We com¬ 
pare the sum of the stellar masses of the cluster member 
galaxies and their total (halo) mass for these 13 clusters, 
and these two mass measures — halo mass (M 200 ) and 
total stellar mass (EM*) — show a positive correlation 
as shown in the left pane of Figure [9j In the middle 
panel of this figure, we show the redshift versus the mass 
ratio (EM*/M 2 oo) of clusters. Here, we cannot see any 
clear redshift-dependent trend, and the mean mass ratio 
is 0.013, shown as a blue horizontal line in this panel. 
The right panel shows the redshift versus M 20 o of our 46 
clusters. Here, red larg e circles are the 13 cluste rs with 
their M 20 o taken from Finoguenov et al.1 ( 201 0). Blue 
small circles are the remaining 33 clusters, for which we 
assign a halo mass assuming the stellar-to-halo mass ra¬ 
tio of 0.013. These clusters have halo masses of 3 x 10 13 
to 2 x 10 14 Mq, and total stellar masses of 4 x 10 11 to 
2 x 10 12 M q . 



Figure 9. Left: Correlation between the cluster halo mass 
(M 200 ) and the total stell ar mass (SM*) of 13 g alaxy clusters 
which are also id entified by IFin ogueno v et al.l (120101 ). Cluster halo 
masses are from [Finoguenov et all 11201 a) . As can bee seen here, 
there exists positive correlation between the total stellar mass and 
the halo mass. Blue line shows the constant £M*/M 200 ratio of 
0.013. Middle: The stellar-to-halo mass ratios and their redshifts 
for the same 13 galaxy clusters. Right: Cluster halo masses and 
their redshifts for all 46 galaxy clusters. Red circles are the same 
13 galaxy clusters shown in left two panels. Blue small circles are 
the remaining 33 clusters. For these 33 clusters, we assume the 
stellar-to-halo mass ratio of 0.013. 


5. COLOR AND STAR FORMATION PROPERTIES 

In this section, we take a close look at the color and SF 
properties of galaxies belonging to the cluster candidates, 
compare these properties with those of field galaxies in 
the same UDS field, and investigate their evolutionary 
trend. 


5.1. Optical Color and its Evolution 
5.1.1. Color of Cluster and Field Galaxies 

Figure ITUl shows (u — g )o color versus stellar mass dia¬ 
grams of member galaxies of each cluster candidate. As 
shown in this figure, the red sequence is formed through¬ 
out the whole redshift range, while the fraction of the 
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Figure 10. The (u — g )o versus M* plots of cluster member galax¬ 
ies. Each panel shows the color-mass diagram for each individual 
cluster candidates. In each panel, the redshift of the cluster is 
given. Red dot-dashed line in each panel shows the dividing line 
((u — g )o = 1.2) between the red and the blue galaxies. 



log M, log M, log M, log M, 

Figure 10. Continued. Black vertical dotted line in each panel 
shows the stellar mass limit (shown as red curves in Figure 5) at 
corresponding redshift. We do not show this limit when it is smaller 
than 9.1. 

galaxies with red color decreases with increasing redshift. 
Hereafter, we refer the galaxies with {u — g)o > 1.2 as red 
galaxies, and the ones with ( u — g) o < 1.2 as blue galax¬ 
ies. In clusters, a significant sequence of red galaxies 
starts to appear at redshift as high as z ~ 1.6. 

At the very highest redshift (z > 1.9), only several 
massive galaxies have red (> 1.2) colors while progres¬ 
sively lower mass galaxies join the red sequence with de¬ 
creasing redshift. At z < 1, we can see that the red se¬ 
quence is well extended down to very low stellar mass (log 
(M*/M 0 ) < 9.5) in most of the clusters. This is a cluster 
version of downsizing phenomenon, which is also shown 
in the study of two gala xy clusters (at z ~ 0.8 and ~ 1.2) 
by iNantais et al.l ( 2013 ). We confirm this trend of mass- 
dependent timing of red sequence formation with a signif¬ 
icantly larger sample. A similar, luminosity-dependent 
tren d in a defic it of r ed galaxies was suggested earlier by 
iGilbank et al.l (l2008h . Please note that we can detect red 
galaxies down to a low-mass limit (= 10 91 Mq) through¬ 
out the redshift range (0.5 < z < 2.0), thanks to the very 
deep optical data from SXDS. Therefore, the dearth of 
low-mass red galaxies at highest redshift bins is real 


i.e., not affected by the magnitude limit of the data. 
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Figure 11. Top: The (u — g)o color distributions of cluster galax¬ 
ies in five redshift bins. The cluster galaxies are summed in each 
redshift bin. Bottom: The (u — ^^-distributions of the field 
galaxies in the same redshift bins as in the top row. The rest- 
frame color is derived from the best-fit BC03 template for each 
galaxy. Red vertical line in each panel shows the dividing line 
((u — g) o = 1.2) between red and blue galaxies. Both in cluster- 
and field-environment, color distributions show clear bimodality, 
except at the highest redshift bin (z ~ 1.85). There is no clear 
difference in color distribution between clusters and field at the 
highest redshift while the distributions are clearly distinguished at 
the lowest redshift (z ~ 0.65). The evolution of color distribution 
in clusters starts to deviate from the field ones from z ~ 1.25. 

In Figure QT] we present the distributions of (u — g )o 
color at five redshift bins (0.5 < 2 < 0.8, 0.8 < z < 1.1, 
1.1 < z < 1.4, 1.4 < z < 1.7, and 1.7 < z < 2.0) for 
cluster (top row) and field galaxies (bottom row). As 
can be seen in the color-M* diagrams, {u — g) 0 color 
shows a clear bimodal distribution, which suggests that 
the color transition of galaxies occurs in a short time 
scale. At redshifts greater than z ~ 1.4, the clusters are 
dominated by blue galaxies, while at z < 0.8, clusters are 
red-dominated — i.e., majority of the cluster galaxies 
are red at z < 0.8. But, red galaxies exist even at the 
highest redshift bin (z ~ 1.85). We can also observe 
the reddening of the peak color of blue galaxies with 
decreasing redshift, which reflects the evolution of the 
average stellar population of blue galaxies with time. 

When we compare the (u — g) 0 color distributions of 
the cluster galaxies (top row) and of the field galaxies 
(bottom row), we can see that the color distribution is 
similar in clusters and in field at the highest redshift 
bin (z > 1.7). At z < 1.4, the color distribution in clus¬ 
ters becomes distinguished from that of the field galaxies. 
Clusters become to be more red-dominated compared to 
field at the same redshift. At the lowest redshift bin, blue 
galaxies are still the major population in field, while red 
galaxies clearly dominate in clusters. From this, we can 
conclude that the color transition of galaxies are acceler¬ 
ated in clusters at redshift lower than z ~ 1.4. 

5.1.2. Color Evolution in Clusters and Field 

In Figure IT21 we show the redshift evolution of the me¬ 
dian color of blue galaxies, both for clusters and field. 
In this figure, the small green circles show the median 
color of blue galaxies and the redshift for each cluster 
candidate, and the red, large circles show the median 
colors and the median redshifts of blue cluster galaxies 
summed in each redshift bin. The blue diamonds rep¬ 
resent the median values of field blue galaxies at each 
redshift bin. The corresponding SDSS value, which is 
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derived from the MPA-JHU catalog^, is shown as the 
magenta circle. The expected color evolution of galax¬ 
ies from the BC03 model with various (delayed) SFHs 
are also shown as the green dotted lines (r = 1.5 Gyr, 
Zf = 5), blue dashed lines (r = 2.0 Gyr, Zf = 5), and red 
solid lines (r = 3.0 Gyr, Zf = 8). Here, 2 / is the forma¬ 
tion redshift. For each color, the lines from the bottom 
to the top are with the increasing amount of dust at¬ 
tenuation from E(B — V) = 0 with the increment with 
0 . 1 . 



Figure 12. The evolution of the median (u — g) o color of blue 
galaxies as a function of the redshift. The small green circles show 
the median colors of the blue galaxies for each cluster candidate. 
The large red circles and the blue diamonds are the median (u — g)o 
values of blue galaxies in cluster (red circle) and in field (blue di¬ 
amond). Vertical error bars show the first and the third quartiles, 
and red horizontal error bars show the range of redshift of individ¬ 
ual clusters included at each redshift bin. The magenta circle shows 
the SDSS value. The lines show the predicted evolution of (it — < 7)0 
colors from BC03 with varying values of r and E(B — V). The 
green dotted lines represent the color evolution with r = 1.5 Gyr 
and Zf = 5, where Zf is the formation redshift. The blue dashed 
lines are for the color evolution with r — 2.0 Gyr and Zf = 5, 
and the red solid lines are for r = 3.0 Gyr and Zf = 8. When 
Zf = 5 (blue and green lines), the age ( t ) is 3.0 Gyr and 4.6 Gyr 
at 2 = 1.5 and z = 1. If Zf =8 (red lines), the ages are 3.6 and 
5.1 Gyr at z = 1.5 and 1.0. Three lines with same color and style 
represent the color evolution of templates with E(B — V) = 0.0, 
0.1, and 0.2 from bottom to top. Solar metallicity is assumed in 
the case of the synthetic color evolution. This reddening of blue 
galaxy color with decreasing redshift indicates the average ageing 
of stellar populations in blue galaxies. 

As already mentioned, the blue peak evolves to be¬ 
come redder with decreasing redshift, reflecting the av¬ 
erage aging of the stellar population in the blue galaxies. 
Apparently, the BC03 models with r = 3.0 Gyr, Zf = 8 
(with E(B — V) = 0.0 and 0.1 - i.e., upper two red 

solid lines) or the models with r = 2.0 Gyr, Zf = 5 
(blue dashed curves) seem to well bracket the observa¬ 
tional trend of blue galaxies. This does not necessarily 
mean that the individual galaxies would evolve with this 
SFH, even though it can be a representation of global 
SFH within this redshift range. Interestingly, these SFHs 
have a peak of SFR at 2 ~ 1.8 or 2, showing a coinci¬ 
dence with the peak of the global SFR density evolution. 
It should be noted that the difference between cluster- 
and field-galaxies is not significant, considering the as¬ 
sociated ranges of error bars and the fact that median 

8 http://www.mpa-garching.mpg.de/SDSS/ 


(u — g) o color can be affected by non-negligible fraction 
of green galaxies. 

In interpreting this color evolution of blue galaxies, we 
should consider the fact that the blue galaxy population 
is not a closed one. As galaxies evolve to become redder, 
these galaxies would move out to the red galaxy popu¬ 
lation. Also, low mass blue galaxies that were originally 
out of our mass-cut would be progressively included into 
the blue galaxy population as their stellar mass grows. 
The effect of these is to prevent the median color of 
the blue galaxy population from being reddened quickly. 
Therefore, it is very likely that the SFRs of the blue 
galaxies would evolve more rapidly (i.e., having smaller 
value of r) than the one represented by the red curves in 
Figure [12l Also, we cannot rule out the possibility of the 
change in the median dust attenuation values of the blue 
galaxies with redshift. The distributions of E(B — V ) 
show little, if any, evolution within the redshift range, 
0-5 < 2 < 2.0, with a hint of slight decrease of the mean 
E(B — V) at the highest redshift bins (z > 1.6). If we ac¬ 
cept this small amount of evolution of dust attenuation, 
it means that the higher dust attenuation at lower red¬ 
shift contributes, to some extent, to the redder color (of 
the blue galaxies) at lower redshift, requiring larger value 
of r. This will compensate, to some extent, the effects of 
migration of fading galaxies out of the blue population, 
making our estimation of the (global representative) SFH 
more robust. 

The evolution of the (u — g) o color of quiescent galaxies 
is shown in Figure [13] Unlike blue galaxies, the quies¬ 
cent galaxy population show little color evolution within 
the redshift range of 0.5 < 2 < 2.0 either in clusters or 
in field. This no-evolution of (u — g) o color reflects the 
fact that new quiescent galaxies (with the colors bluer 
than the existing quiescent ones) keep joining the quies¬ 
cent population while the color of the existing quiescent 
galaxies becomes redder with time. The dashed curves 
with various colors in Figure [13] are the evolutionary path 
of BC03 quiescent galaxies with various formation red¬ 
shift. These curves are shown to guide how the (u — g) o 
colors would evolve if no additional quiescent galaxies 
with bluer color join the existing quiescent population. 
Purple circles are the corresponding SDSS colors in var¬ 
ious stellar mass bins - log(M*/M Q ) = 9.7 for the bot¬ 
tom one to log(M*/M q) = 11.1 for the top one with 
an increment of 0.2 in logarithmic scale. Analysing red 
galaxy spectra at z ~ 0.9 from t he Deep Extragalactic 
Evolutionary Probe 2 (DEEP2), iScliiavon et all ()2006l ) 
found relatively young ages (~ 1 Gyr) for these red 
galaxies with little or no emission lines. One of their 
interpretation is that blue SF galaxi es keep j oining: this 
red population co ntinuously (see also lHarker et al .1120061 
iFaber et, al.l l2007h maintaining their mean stellar popu¬ 
lation ages young. 

5.1.3. Red Star-forming Galaxies 

The (u — g)o color and sSFR of galaxies show a broad 
correlation in a sense that redder galaxies have in general 
lower sSFR. However, there are a non-negligible fraction 
of galaxies whose (u — g) o color is red (> 1.2) but are 
still forming stars actively. As shown in Figure [TT! these 
red star-forming (SF) galaxies are, on average, dustier, 
older, and have lower sSFR than blue SF galaxies. No 
red SF galaxy has sSFR > 10 -8 yr -1 to be classified as 
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Figure 13. The evolution of the median (u—g )o color of quiescent 
galaxies as a function of redshift. The symbol assignment is same 
as Figure IT2] (u — g )o color of quiescent galaxies shows little evo¬ 
lution within the redshift range, 0.5 < z < 2.0. This implies that 
newly quenched galaxies keep joining the quiescent galaxy popu¬ 
lation within this redshift range. Dashed curves are the evolution¬ 
ary tracks of quiescent galaxy templates from BC03 with various 
formation redshift ( Zf ) as indicated in the figure. Purple circles 
are the corresponding colors of SDSS galaxies with different stellar 
masses, which increase from bottom (log(M*/MQ) ^ 9.7) to top 
(log(M* /Mq) ~ 11.1) with a decrement of 0.2. Bracketed numbers 
are the median stellar masses at corresponding redshifts for field 
and cluster galaxies, respectively. 

‘starbursts’. These indicate that red SF galaxies are in 
the fading stage, and migrating into the red quiescent 
population. 
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Figure 14. Distributions of dust extinction (E(B-V); Left), mean 
stellar population age (Middle), and sSFR (Right) estimated 
from SED-fitting for blue SF (blue hatched histogram) and red 
SF (red filled histogram) galaxies. On average, red SF galaxies are 
more dustier as well as older with lower sSFRs. 


We investigate the morphological properties of red 
SF galaxies and compare them with the other galaxy 
populations to see how many red SF galaxies could be 
starburst galaxies resulting from gas-rich (wet) merger. 
A portion of the UDS region (~ 22.3' x 9') is cov¬ 
ered by the CANDELS (Cosmic Assembl y Near-i n frared 
Deep Extragalactic Le gacy Survey; iGrogin et al.1 120111 : 
iKoekemoer et, alll2011fl program, and observed with the 
WFC3 (Wide Field Camera 3) on board the HST. From 
the deep WFC3 H-band (F160W) image, which reaches 
to the 5-cr magnitude limit of 27.45 with FWHM of 0.2" 
(jGalametz et al.l 1201 3). the Gini coefficients, G, as well 
as the M 2Q parameters of the galaxies are measured. Fig- 
ureHHshows the distributions of the Gini coefficients and 
the M 2 o values of the blue, the red SF, and the red qui¬ 
escent galaxies. The dividing lines betwee n the early -, 
late-type and major mergers are from lLotz et al.l (|2008f ). 
There are only a few merger candidates among red SF 


galaxies (in the middle panels, represented as green sym¬ 
bols) and most of the red SF galaxies lie around the 
boundary between blue (left panels) and red quiescent 
galaxies (right panels), supporting the idea that these 
red SF galaxies are in transition phase from blue galax¬ 
ies to red quiescent ones. 
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Figure 15. The G-M 20 diagrams of galaxies in the sub-region of 
the UDS covered by the CANDELS. Left column (blue points): 
Blue galaxies. Middle column (green points): Red SF galaxies. 
Right column (red points): Re d quiescent galaxies. The dividing 
lines are from [Lotz et al.l ll200811 . Majority of red SF galaxies lie 
in the region between blue galaxies and red quiescent ones in this 
diagram, and only a few red SF galaxies are merger candidates. 


We also visually inspect the HST/ WFC3 images of 
our galaxies. The blue SF population mostly consists of 
extended objects with disk-like morphology with small 
fraction of disturbed ones hinting the recent experience of 
merger or interaction. On the other hand, red SF galax¬ 
ies are mixture of disky, spherical, and compact objects, 
reinforcing that these objects are in transition from blue 
galaxies into red quiescent ones. In Figure 1161 we show 
example HST /WFC3 17-band images of cluster galaxies 
at z ~ 0.65. 

5.2. SF properties and Evolution of Cluster and Field 

Galaxies 

In previous section, we investigated the rest-frame op¬ 
tical color properties of galaxies in clusters and field. 
Now, in this section, we concentrate on the SF properties 
of galaxies and analyse how the SF properties of galax¬ 
ies within the high-redshift clusters evolve with time and 
also compare this with that of the field galaxies in the 
same redshift range. 

In Figure fTTl we compare the distributions of sSFR of 
the cluster- and the field galaxies at five redshift bins 
with the same redshift binning as in the Figure 1111 In 
this figure, we assign sSFR = 10^ 12 yr -1 to all quies¬ 
cent galaxies with sSFR < 10 -12 yr -1 . Similarly with 
color distribution, the sSFR distribution of SF galaxies 
in clusters and in the field are similar at high redshift 
bins (z > 1.4), while it shows a clear difference from that 
of the field galaxies at the lowest redshift bin (z < 0.8). 
For the quiescent galaxy fraction (over SF ones), the dif¬ 
ference between clusters and field starts to show up from 
2 ~ 1.25. While the evolutionary trends of sSFR dis¬ 
tributions is similar with that of the color distributions 
shown in Figure [Tl] the sSFR distribution changes more 
slowly (especially for the cluster galaxies) than the color 
distribution. For example, at 2 ~ 0.95, the color distri¬ 
butions are already different between cluster and field, 
while the sSFR distributions of SF galaxies still remain 
similar in field- and cluster-environments. 
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Figure 16. HST/ WFC3 R-band (F160W) images of cluster galaxies at z ~ 0.65. (Left) Examples of blue star-forming galaxies. 
(Middle) Red star-forming galaxies. (Right) Red quiescent galaxies. Each image cut size is ~ b" X 5". We show the values of Gini 
coefficient (G) and M20 of each galaxy as well. 
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Figure 17. Top: The sSFR distributions of cluster galaxies in 
five redshift bins. The cluster galaxies are summed in each redshift 
bin. Bottom: The sSFR distributions of the field galaxies in the 
same redshift bins as in the top row. The sSFR is given in the 
unit of Gyr - b In each panel, we assign sSFR = 10 —3 Gyr -1 to all 
quiescent galaxies with sSFR j ICR 3 Gyr -1 . The sSFR, distribution 
of cluster SF galaxies shows clear difference from that of field SF 
galaxies only at the lowest redshift bin (z ~ 0.65), while relative 
fraction of quiescent galaxies shows difference between clusters and 
field at 2 < 1.4. 


Next, we compare the evolution of the quiescent galaxy 
fraction in cluster and in field in Figure IT51 dividing the 
sample into three stellar-mass bins : 9.5 < log M*/Mq 
< 10.0 (left panel), 10.0 < log M*/Mq < 10.5 (middle), 
and log M*/M 0 > 10.5 (right). In this figure, the filled 
blue diamonds represent the quiescent fraction of field 
galaxies at each redshift bin, while the quiescent frac¬ 
tion of individual clusters is shown as magenta dots. At 
several redshift bins, we sum the cluster galaxies in each 
2 -bin, and calculate the mean and the standard devia¬ 
tion. This is shown as the red diamonds with error bar 
in the figure. The red and blue circles are for the S PSS 
(Sloan Digital Sky Survey) galaxies from iBaldrv et alJ 
( 2006 ). We derive these values from their Equation (9), 
and choose the values with a = —0.3 as field values and 
a = 0.9 as cluster values in that equation. An d, we ap ply 
the co rrection for red SF galaxies based on lHaines et alJ 

(l2008h . 

From Figure 1181 we find several interesting aspects in 
the evolution of quiescent galaxy fraction. First, the in¬ 
crease of the quiescent fraction slows down from the red¬ 
shift 2 ~ 1.3 —1.4 for galaxies with their stellar mass, log 
(M*/M 0 ) > 10, more significantly for the field galax¬ 
ies, while it evolves fast from z ~ 2 down to z ~ 1.3 
for massive galaxies (M* > 1O 1O ' 5 M 0 ). At the two high 
stellar-mass bins, the quiescent fraction remains nearly 
unchanged from z ~ 1.3 — 1.4 both in held and in clus¬ 
ter. Comparison with the local values (the red and the 
blue circles) shows that the quiescent fraction has already 
reached to the local value for the most massive galax¬ 
ies (log (M*/M 0 ) > 10.5) at z ~ 0.6 both in clusters 
and in held (the right panel), while the growth of qui¬ 
escent fraction must be accelerated at 2 < 0.5 to match 
the SDSS values in the case of low mass galaxies with 
log (M*/Mq) < 10. This near-constant quiescent frac¬ 
tion of massive galax ies is in good agreement with th e 
previous studies (e.g.. llm et a 1.1 120021 : Ilbert et aI1l2013D . 
This change in the increase of quiescent fraction indi¬ 
cates that the redshift range 2 > 1.3 dehnes the era of 
the rapid build-up of massive quiescent galaxies, which 
is in broad agreement with the hnding of actively star¬ 
forming galaxy clusters at 2 > 1.3 by several authors 
(e.g., Tran et al. 2010; Zeimann et al. 2012; Santos et 
al. 2014). 




Redshift Redshift Redshift 

Figure 18. Upper: The quiescent fraction evolution of galax¬ 
ies in three stellar-mass bins (log M*/Mq ~ 9.75 (left), 10.25 
(middle), and > 10.5 (right). The purple dots show the quies¬ 
cent galaxy fraction of the individual cluster candidate and the 
red diamonds with the error bar are the number-weighted mean 
and the standard deviation of cluster galaxies summed at discrete 
redshift bins (A z = 0.3). The blue diamonds are the quiescent frac¬ 
tion of the field galaxies at each redshift bin (Az = 0.2). The red 
and blue circles in each panel show the cluster and field quiescent 
fract ion at each correspo nding stellar mass of the SDSS galaxies 
from [Baldrv et ahl Q2006A . The difference of quiescent fraction in 
different stellar mass bins is greater than the difference between 
different environment in the same stellar mass bin. Solid curves in 
each panel show the best-fit quiescent fraction evolution (Equation 
4) for cluster (red) and field (blue) galaxies, respectively. Lower: 
The excess of the quiescent fraction in the clusters over the field 
value, defined as / q,cluser /fq,fieldi where fq,cluster &nd fq,field 
are quiescent fraction in the clusters and in the field, respectively. 
The thick error bars show the standard deviation at each redshift 
bin, while the thin error bars reflect the spread among clusters in 
the given redshift bin. We observe this excess becomes significant 
(> 1.5) only in the smallest mass bin (lower left panel) at z < 1.2. 


Second, the fraction of quiescent galaxies is higher 
in clusters than in field at 2 < 1.4 for galaxies with 
M* < 10 1o ' 5 Mq, even though the scatter among the in¬ 
dividual cluster candidates is quite large. From this, we 
can speculate that the environmental quenching becomes 
to work more strongly at redshift, 2 < 1.4, but with a cer¬ 
tain amount of cluster-to-cluster variation. Also, we can 
see that the difference in the quiescent galaxy fraction 
between the cluster and the field environments is more 
significant for low-mass (log (M*/M 0 ) < 10) galaxies 
(the left panel) at redshift 2 < 1. 

Last, as can be seen clearly in this figure, the quiescent 
fraction is a strong function of the stellar mass rather 
than their environment. At the highest mass bin, the 
cross-over, which we define as the redshift or the epoch 
when the quiescent fraction starts to exceed 50 %, occurs 
already at 2 ~ 1.5. In sharp contrast, in the lowest mass 
bin, the quiescent fraction never reaches the cross-over 
down to 2 ~ 0.5 both in field and in cluster. Actually, 
the quiescent fraction is < 0.2 throughout the redshift 
range in field. This strong mass-dependence of quiescent 
fraction evolutio n well agrees with the recent results of 
iMoustakas et al.l (120131) — who found a nearly constant 
number density of massive quiescent galaxies at 2 < 1 
while the corresponding value rises rapidly with decreas¬ 
ing redshift for less massive galaxies, and of lHuang et all 
(120131 ). This strong mass dependence and the weaker 
environmental dependence dictates that the stellar mass 
plays a more dominant role in shaping the SFH of galax¬ 
ies than the cluster-specific processes. 
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To gain an insight how rapidly (or how slowly) star 
formation is quenched for galaxies with different stellar 
masses, we devise a simple model explaining the evolu¬ 
tion of star-forming galaxy fraction, f s f(t ), as follows. 

fsf(t) = fsf( O)x(i)^, (3) 

where, f s f( 0) is the fraction of SF galaxies at an initial 
time t = 0. In this model, 7 is a half-life of SF galaxy 
population, which means that the fraction of SF galaxies 
becomes half of its initial values after the time t = 7 . 
This simplified model assumes that the total number of 
galaxies does not change with time - i.e., no galaxy is 
added to or removed from the sample, and no merger 
occurs between galaxies. 

Then, the quiescent fraction, f qu i es . would increase as, 

W(t) = l-/a/(0)x((4) 

We fit our data points to this simple model, and we 
show this fquies. for the best-fit values of 7 ’s as solid 
curves in each panel in the upper row of Figure 1181 In 
each panel, the red and the blue curves are for cluster and 
field galaxies, respectively. For cluster galaxies at 0.5 < 
z < 2.0 (solid red curves), the best-fit values for 7 are 
15±2.2, 9.7±3.6 and 4.4±1.9 Gyr for the lowest, middle 
and highest stellar-mass bins. In the case of field galaxies 
(solid blue curves), the corresponding values are 46±7.1, 
18±3.9, and 5.2±1.4 Gyr, each. This again demonstrates 
that: ( 1 ) massive galaxies become quiescent more rapidly 
both in clusters and in field (i.e., having smaller value of 
7 ) than less massive galaxies, and ( 2 ) star formation is 
more rapidly quenched in clusters than in field. 


0.5 <z<1.0 1.0<z<1.5 1.5 <z<2.0 



Figure 19. The quiescent fraction evolution of galaxies in three 
redshift bins (z ~ 0.75 (left), 1.25 (middle), and 1.75 (right). 
The red and blue diamonds with the error bar are the number- 
weighted mean and the standard deviation of cluster and field 
galaxies summed at discrete stellar-mass bins, respectively. The 
clear dependence of quiescent fraction on the stellar-mass can be 
seen. 

The lower panels of Figure |T 8 ] show the ratio of the 
quiescent fraction in clusters ( f q ,dustesr ) to that in the 
field ( f q ,field, )■ Here, we take an average value of the 
highest three redshift bins (shown as the black circles) 
and show it as a red diamond in each panel. In the two 
high mass bins (middle and right panel), the fraction is 
always < 1.5 with little evolution. However, this fraction 
increases rapidly with decreasing redshift for low mass 
galaxies (log M*/M 0 < 10 ), and is higher than that in 
the higher mass bins at z < 1.4. This also indicates that 
the environmental quenching becomes more significant 
for low mass galaxies at z < 1.4. 

To see the effect of stellar mass in determining the 
quiescent fraction more clearly, we show the stellar-mass 


dependent quiescent fraction at three redshift bins in Fig- 
ure[l9l In this figure, we can clearly see the strong stellar- 
mass dependence of the quiescent galaxy fraction at all 
three redshift bins as well as a clear excess in the quies¬ 
cent galaxy fraction in clusters compared to the field for 
low-mass galaxies (log (M*/M 0 ) < 10.5) at the lowest 
redshift bin (left panel). 

Summarizing, our investigation of the quiescent galaxy 
fraction and its evolution reveals several important as¬ 
pects about quiescent galaxy formation: ( 1 ) At z > 1.3, 
quiescent, massive galaxies were built up rapidly, with 
a similar rate in both clusters and field. By z 1.3, 
most of the massive,quiescent galaxies were built up. 
( 2 ) At z < 1.3, the environmental dependence of qui¬ 
escent fraction become visible. The difference between 
cluster- and field-environment is clearer for low-mass (log 
(M*/M 0 ) < 10.0) galaxies. (3) Stellar mass plays more 
dominant role in determining the quiescent fraction than 
the environment throughout the entire redshift range. 

5.3. Quenching Efficiency 

In the previous section, we have shown that the 
quiescent galaxy fraction shows a clear stellar-mass- 
dependent trend, in a sense that the fraction of quies¬ 
cent galaxies is higher for more massive galaxies through¬ 
out the redshift range, 0.5 < 2 < 2 (upper panels 
of Figure USD, and also that the excess of the quies¬ 
cent galaxy fraction in clusters over field environment 
is higher for low-mass (< 1O 1 O M 0 ) galaxies, especially 
at redshifts lower than z < 1.2 (lower panels of Fig¬ 
ure USD- Now, we investigate the environmen tal quenc h¬ 
ing efficiency (e.g.lvan den Bosch et al.ll2008t IPeng et al.1 
l201?f IQuadri et alJI2012H ~ We define this environmental 

quenching efficiency as, 

( fq, cluster fq, field) / fsf^fieldt ( 5 ) 

where f q duster and f q j ie id are the fraction of quies¬ 
cent galaxies in clusters and in field, respectively, and 
fsf,field is the fraction of star-forming galaxies in field 
(= 1 fq,field)• This quantity measures the fraction of SF 
galaxies in field that would have become quiescent if they 
were in cluster. Figure [20] shows the evolution of the en¬ 
vironmental quenching efficiency in different stellar mass 
bins. Our results shows no clear stellar-mass dependent 
trend - even though the scatter is larger in the highest 
mass bin due to the small number of massive galaxies 
and the small f s f,field v alues in this mass bin — con¬ 
firmi ng previous results (jPeng et al.l 120101 : IQuadri et al.1 
I 2012 T) but extending the probed redshift range to z ~ 2 . 

How can this result — no significant mass dependence 
in environmental quenching efficiency — be reconciled 
with the results in the previous section (e.g., Figure 1151) 
the high excess of the quiescent fraction in clusters in 
the lowest mass bin? This difference arises because there 
is another quenching mechanism — i.e., mass quenching. 
Massive galaxies are affected strongly by mass quench¬ 
ing, and they are mostly quiescent already at z ~ 1.3, 
leaving little room for the environmental quenching to 
play a significant role. On the other hand, the lower 
mass galaxies are quenched only mainly through the en¬ 
vironment effects that becomes significant at 2 < 1.3. 

6. SUMMARY AND CONCLUSION 
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Figure 20. Evolution of environmental quenching efficiency, de¬ 
fined as the excess of quiescent galaxy fraction in clusters over field 
divided by SF galaxy fraction in field. Three panels show this evo¬ 
lution in different stellar mass bins. Red circles are our results and 
blue circles are SDSS values. There is no clear difference between 
different stellar mass bins, except large scatter in the highest mass 
bin (log ( M*/Mq ) > 10.5; right panel). 


In this paper, using deep optical to MIR data in 
the UKIDSS/UDS field, we have found 46 high-redshift 
galaxy clusters up to 2 < 2, among which 27 are newly 
found. We analyse the stellar population properties, such 
as color and SFR, of galaxies both in cluster and field en¬ 
vironments over a wide redshift range to understand the 
effects of cluster environment on the galaxy evolution. 

Through this analysis, we have found that the qui¬ 
escent galaxy fraction increases rapidly with decreasing 
redshift at z > 1.3, and the increase is slowed down at 
lower redshift for massive (log (M*/Mq) > 10.5) galax¬ 
ies. This trend points to the redshift range 2 > 1.3 as 
the era of the rapid build-up of massive quiescent galaxies 
— the epoch when many massive galaxies stop their star 
formation and become quiescent. The difference in the 
quiescent galaxy fraction between clusters and the field 
increases at z < 1.2-1.4, which coincides with the epoch 
when the increase in the quiescent galaxy fraction begins 
to slow down. In this redshift range (0.5 < z < 1.4), 
the quiescent galaxy fraction remains nearly unchanged 
in field, while it keeps increasing in clusters (but more 
slowly than at z > 1.4). This difference is only signifi¬ 
cant for low mass galaxies with log (M*/Mq) < 10.0. 

The environmental quenching efficiency shows no clear 
stellar-mass depen dence, which is in agreement with pre - 
vious results (e.g., iPeng et al.ll2010HQuadri et alJl2012h . 
while our results extend this up to higher redshift (z ~ 
2). However, the effects of environmental quenching ap¬ 
pears more significantly for low-mass galaxies, because 
massive galaxies are affected by another quenching mech¬ 
anism — i.e., mass-quenching — and most of these mas¬ 
sive galaxies are already quenched at z > 1. 

At z > 1.4, the quiescent fraction of galaxies shows no 
clear dependence on their environment, in contrast to its 
strong dependence on stellar mass. This infers that the 
SFH of galaxies is mainly shaped by their stellar mass 
during the early phase of evolution. The effects of the 
cluster environment on determining the quiescent frac¬ 
tion of galaxies seem to be more significant for less mas¬ 
sive galaxies with M* < 10 10 M & , for which mass quench¬ 
ing does not yet affect their SF activity, at z < 1.4. The 
delayed appearance of the cluster-environment effects on 
the quiescent galaxy formation for low-mass galaxies im¬ 
plies that the cluster-specific processes which are mainly 
responsible for the excess of the quiescent fraction in clus¬ 
ters compared to the field region may be gradual, time- 
taking processes — like the strangulation (|Larson et al.1 


1980; iBalogh et al.l I2000H — rather than abrupt or vi¬ 
olent processes - like ram-pressure stripping, for low- 
mass galaxies. The investigation of the morphology of 
red SF galaxies — which shows that the most transition 
from the blue SF galaxies to red quiescent ones in clus¬ 
ters does not affecting their morphological appearance 
significantly — supports this idea. 

One interesting question is why the SFHs of galax¬ 
ies are affected by their stellar mass. More specifically, 
why is the quiescent galaxy fraction higher for more mas¬ 
sive galaxies at all environment over a wide range of 
redshift (z < 2) and without strong environmental de¬ 
pendence? One possible answer is that the dynamical 
time-scale would be shorter for the galaxies with greater 
stellar masses than less massive galaxies. If the SFR of 
galaxies depends on this dynamical time-scale, in a sense 
that shorter dynamical time-scale leads to higher SFR 
- as assumed in man y galaxy formation models (e.g., 
iSomerville et alj|2012l f. more massive galaxies will con¬ 
sume their gas more rapidly to arrive the red quiescent 
galaxy sequence earlier than low mass galaxies. Major 
mergers between massive SF galaxies can also accelerate 
this fast gas consumption. Another possibility is that any 
negative feedback which works preferentially for massive 
galaxies — such as AGN feedback (e.g. jHopkins et al.1 
12006b ISomerville et al .1120081 iSantos et al.l 1201 -'ll) or halo 
mass quenching dBirnboim fc Dek cl 2003) — is the driver 
of the mass-quenching phenomena. 
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Figure Al. The quiescent fraction evolution of galaxies in three 
stellar-mass bins (log M*/Mq ~ 9.75 (left), 10.25 (middle), and 
> 10.5 (right) when we apply constant sSFR cut (sSFR = 10 —10,7 
yr —1 ). The sym bol assignment is same as in Figure [THl Similarly 
with Figure [lH] the difference of quiescent fraction in different stel¬ 
lar mass bins is greater than the difference between different envi¬ 
ronment in the same stellar mass bin in both cases. 


One of the interesting results of this work is that the 
quiescent galaxy fraction drops rapidly as we approach 
the redshift, 2 ~ 2, independent of the environment in 
which they reside. This result indicates that the major¬ 
ity of galaxies are actively forming stars at z> 2. Com¬ 
bined with the fact that the quiescent fraction within 
clusters and in the field are similar at this epoch, this in¬ 
dicates that studies of clusters or proto-clusters at z > 2 
would be crucial in revealing the very initial properties 
of the forming galaxy clusters as well as the evolution of 
galaxies in these massive structures. Also, we can expect 
that any cluster finding methods using the presence of 
old and quiescent galaxies within clusters (for example, 
like red-sequence technique) would miss many clusters 
or proto-clusters at z > 2. In our work, both of the 
identification of the galaxy clusters at z > 1.5, as well 
as the reliable estimation of the SFR of galaxies through 
the SED-htting are possible thanks to the deep NIR data 
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from the UKIDSS. This implies that reliable NIR (photo¬ 
metric or spectroscopic) data are essential for the future 
search and the study of the high-redsliift galaxy clusters 
during this important epoch (z > 1.5). Therefore, we can 
expect a big leap in our understanding of the properties 
and the evolution of the high-redsliift galaxy clusters or 
proto-clusters as well as of the galaxy evolution with the 
near-future NIR facilities either in space — like JWST 
(James Webb Space Telescope) — or on the ground — 
like GMT (Giant Magellan Telescope) with the GMACS 
(with the NIR-arm) and later with the NIRMOS. 
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This work was supported by the National Research 
Foundation of Korea (NRF) grant, No. 2008-0060544, 
funded by the Korea go vernment (MSIP ). The UKIDSS 
project is defined in lLawrence et all (12007H . The 
UKIDSS uses the UKIRT Wide Field Camera (WFCAM; 
ICasali et al.112007 ). Th e photometric system is described 
in lHewett et al.l (2006I1. and the calibration is described 
in iHodekin et al. (I2009f) . The pipeline processing and 
scie nce a rchi ve are descri bed in Irwin et al (2009, in prep) 
and lHamblv et all ( 2008 ). Part of this work is based on 
observations taken by the CANDELS Multi-Cycle Trea¬ 
sury Program with the NASA/ESA HST, which is op¬ 
erated by the Association of Universities for Research in 
Astronomy, Inc., under NASA contract NAS5-26555. 

Facilities: UKIRT (WFCAM), Subaru (SUPRIME- 
CAM), HST (WFC3). 


Figure A2. Evolution of environmental quenching efficiency, de¬ 
fined as the excess of quiescent galaxy fraction in clusters over 
field divided by SF galaxy fraction in field with constant sSFR cut 
(sS FR = 10 —10,7 yr — - 1 ). The symbol assignment is same as in Fig¬ 
ure [20] Three panels show this evolution in different stellar mass 
bins, and there is no clear difference between different stellar m^^PENDIX 

bmS ' EFFECTS OF DIFFERENT SSFR CUT 

Since galaxies selected with our redshift-dependent cut of sSFR < 1/3 1 are not quite “quiescent” (i.e., no SF) at 
high redshift, we examine here how a constant sSFR cut changes our main conclusion regarding SF quenching. Here, 
we choose a cut of sS FR = 10~ 107 yr -1 , which corresponds to the cut for local quiescent galaxies (e.g. iGallazzi et ahl 
l2n09tlKo et. ahllml . 

Applying this constant sSFR cut, we first analyse the evolution of quiescent galaxies fraction. In Figure IA1I 
which corresponds to upper panels of Figure 18, we show the evolution of quiescent galaxy fraction in different mass 
bins. Not surprisingly, the quiescent galaxy fraction is overall lower than when we apply the l/[3t(z)]-cut. The most 
significant difference is that the quiescent galaxy fraction shows steeper increase from z ~ 2 to z ~ 1.3, for galaxies 
with M* > 10 10 Mq (middle and right panels), especially for M* > 10 10 5 Mq, than in Figure 18. This indicates that SF 
quenching (to local q uiescen t gal axy level) is faster for more massive galaxies. This fast increase of quiescent fraction 
at high-z agrees with iDominguez Sanchez et al.1 (120111) result. 

While there are interesting differences in quiescent fraction evolution with different sSFR cuts, one important trend 
seems to hold: the quiescent galaxy fraction s hows clearer dependence on their stellar mass than the environment. 
The environmental quenching efficiency (Figure[A2j is not mass-dependent, but rises toward lower redshift (at z < 1), 
similarly with Figure 20. 
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Table 1 

Candidate high-redshift galaxy clusters in the UDS 


RA 

(1) 

dec. 

(2) 

z 

(3) 

N gal 

(4) 

T.M * 
(5) 

°OD 

(6) 

Reference 

(7) 

34.48571 

-4.88290 

0.506 

46 

7.08 

7.81 

1 

34.66831 

-5.05773 

0.511 

39 

10.42 

7.25 


34.05109 

-4.74287 

0.534 

42 

5.61 

6.30 


34.70426 

-5.14720 

0.538 

34 

6.80 

5.61 


34.54121 

-5.36501 

0.607 

55 

6.44 

4.81 

2 

34.54105 

-5.26098 

0.626 

63 

10.48 

6.37 


34.35041 

-5.41133 

0.629 

57 

6.88 

5.74 

1 

34.19640 

-5.15057 

0.633 

59 

7.57 

4.89 

1,2 

34.39710 

-5.22284 

0.639 

93 

16.34 

13.60 

1,2 

34.60152 

-5.41210 

0.644 

79 

14.49 

8.78 

1 

34.51734 

-5.52098 

0.648 

46 

6.90 

6.50 


34.63576 

-4.96694 

0.659 

39 

5.20 

5.23 

2 

34.47932 

-5.45278 

0.677 

60 

11.63 

10.05 

1 

34.74237 

-5.12341 

0.684 

32 

6.54 

6.37 


34.37319 

-4.68903 

0.685 

29 

3.81 

6.75 


34.84291 

-4.82145 

0.744 

28 

3.31 

5.49 


34.45845 

-5.50699 

0.774 

42 

5.43 

6.73 


34.42685 

-5.09312 

0.791 

37 

7.44 

5.60 


34.52339 

-4.73828 

0.839 

44 

9.04 

5.73 

1 

34.82770 

-5.08506 

0.865 

39 

9.12 

5.92 

1 

34.63694 

-5.01183 

0.869 

72 

13.25 

9.45 

1 

34.16359 

-4.73395 

0.889 

46 

7.78 

5.12 


34.84166 

-4.88236 

0.890 

53 

8.08 

5.12 

1 

34.34870 

-5.20672 

0.909 

54 

9.74 

7.36 


34.05085 

-4.87705 

0.931 

54 

8.41 

6.16 


34.06289 

-4.71558 

0.932 

51 

6.93 

6.50 


34.53990 

-5.01259 

0.935 

56 

7.96 

5.87 

2 

34.03834 

-5.10964 

1.042 

43 

8.61 

5.18 


34.29435 

-4.79312 

1.042 

27 

6.09 

5.09 

2 

34.59195 

-4.97802 

1.038 

38 

8.43 

5.93 


34.52635 

-5.01766 

1.064 

38 

10.51 

5.68 

2 

34.28226 

-4.82255 

1.109 

34 

9.19 

5.68 

1 

34.61517 

-4.69712 

1.214 

35 

8.57 

4.89 


34.69049 

-4.71229 

1.238 

46 

10.41 

5.34 


34.80897 

-4.93465 

1.246 

43 

6.03 

5.73 


34.53326 

-5.01095 

1.281 

41 

12.55 

6.12 

3 

34.85731 

-4.86560 

1.364 

49 

15.61 

7.69 


34.84134 

-4.72292 

1.366 

32 

8.73 

5.63 


34.53797 

-5.01744 

1.500 

36 

9.48 

5.60 

1,3 

34.07070 

-5.00289 

1.517 

36 

8.18 

5.49 


34.81516 

-4.74594 

1.538 

41 

7.15 

6.18 


34.59165 

-5.16940 

1.640 

66 

14.37 

6.93 

1,4 

34.17832 

-5.15493 

1.664 

62 

11.41 

5.34 


34.80031 

-4.73156 

1.671 

77 

13.34 

6.17 


34.73531 

-5.04661 

1.916 

40 

4.73 

6.89 


34.72378 

-5.17741 

1.983 

35 

6.79 

5.30 



Note. — 

(1) RA in degree 

(2) Declination in degree 

(3) Redshift 

(4) Number of galaxies within 1 Mpc radius from the cluster 
center 

(5) Sum of the stellar masses of the member galaxies in le+11 

Mq 

(6) Overdensity measure, ( N — N)/<jn 

(7) Reference list: 

1. Finoguenov et al. 2010, MNRAS , 403, 2063 

2. van Breukelen et al. 2006, MNRAS , 373, L26 

3. van Breukelen et al. 2007, MNRAS , 382, 971 

4. Papovich et al. 2010, ApJ , 716, 1503 



